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Abstract. At zero temperature, the classical antiferromagnetic Ising model on 
the pyrochlore lattice is a spin disorder phase of the critical spin correlation. It is a 
deconfined phase in that the binding energy of the monopole-anti-monopole pair is 
independent of their distance of separation. We show that turning on a transverse 
magnetic field turns it into the cooperative paramagnet, and the spin correlation 
becomes exponential decay. Furthermore, it introduces the quantum confinement 
(of magnetic monopoles) , where the binding energy of the pair is proportional to 
their distance of separation. This disorder state undergoes adiabatic transition to 
the paramagnetic state in the large field limit. The cff'cctivc Hamiltonian (without 
magnetic monopoles) in small field is the Ising Hamiltonian plus ring exchange 
interax;tion. 
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Quantum number fractionalization is one of the most intriguing topics in physics. 
In particular, searching for the magnetic monopole has been of persistent interest for 
centuries. This interest was recently aroused again by the discovery of the magnetic- 
monopole-like spin excitations in the spin ice systems [2 [5]. The energy of the 
monopole-anti-monopole pair was shown to be inversely proportional to their distance 
of separation, i.e. the Coulomb phase of the magnetic monopole [3]. 

The spin ice systems are found in the pyrochlore lattice. It is a 3-dimensional 
system in which tetrahedra are stacked by sharing their corners as shown in Fig. ([2^) . 
The ground state of the spin ice system is that the spin orientations in each tetrahedron 
are pointing " 2-in-2-out" of the tetrahedron [H [3] , similar to the ice rule [S] . For a 
given classical spin configuration in the ground state manifold, one can create a spin 
excitation by flipping a spin. It results in breaking the ice rule at two neighboring 
tetrahedra. In one of them, the spin orientations will be "l-in-3-out" , giving rise 
to a monopole, and in the other they will be "3-in-l-out" , giving rise to an anti- 
monopole. One can try to hop the anti-monopole by flipping another spin pointing 
inward to be outward in the anti-monopole tetrahedron. If only the Ising interaction 
is taken into account, the distance of separation of the monopole and anti-monopole 
can increase with no cost in energy. In other words, the pair is deconfined because the 
binding energy does not depend on their distance of separation. When considering 
the classical dipolar interaction additionally, the binding energy becomes inversely 
proportional to their distance of separation [3]. 

Unlike ordinary electrons, the emergent magnetic monopole and anti-monopole 
in the spin ice systems are classical particles [3j. It is important to investigate 
whether their quantum counterparts are deconfined or not. Searching for the quantum 
deconfinement and quantum spin disorder phase, namely the spin liquid, has attracted 
tremendous attention in the last few decades [HI HI [S [13 EI]- The questions 
are, however, extremely difficult in the 3-dimensional system, because the quantum 
fluctuation is largely suppressed in three dimensions and, in addition, the limits on 
computational power also restricts the calculation from being conclusive. In this 
Letter, we study the antiferromagnetic quantum Ising model 



where (ij) sums over all nearest neighbor spins, i sums over all spins, J > 0, and S'^ 
are the Pauli spin matrices for spin-1/2. We will numerically demonstrate that the 
ground state of Eq. ([!]) is spin disorder with an exponential decay correlation, i.e. 
disorder from disorder [T^ [T3| . The binding energy of the monopole-anti-monopole 
pair will be calculated by using the generalized random phase approximation. In 
the small field limit, it is ~ '"^j ^■'^d ~ I K in the large field limit, where 
I is the number of bonds separating the pair. In other words, due to a non- 
vanishing transverse field, the magnetic monopoles undergo an abrupt transition 
into the quantum confinement regime. Moreover, we will argue that the new spin 
disorder phase and the paramagnetic phase in the large field limit can be connected 
adiabatically by a continuous transition, and therefore the spin disorder phase at small 
field is a cooperative paramagnet. In conclusion, our work serves as the simplest model 
to host the spin disorder state in the 3-dimensional systems. Finally, the effective 
Hamiltonian in small field will be derived for the connection of our work to the recent 
progress in frustrated spin systems. 

We show the disorder nature by computing the spin correlation length of T{i) = 




(1) 
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Figure 1. (Color online) (a) The spin correlation length i;{L, T, K) as a function 
of 1/L for three sets of T and K. Both T and K are in units of J, and 
the correlation length is in units of the bond length. The spin correlation is 
exponentially decayed and the correlation length is around one bond distance 
down to the temperature as low as 2.5% of K for K = 0.05J and 5% of K 
for K = 0.025J and 0.005J. The weak dimensional dependence indicates that 
our calculations are sufficient for the thermodynamic limit, (b) The ^(T,K) as a 
function of T for three different K values. The correlation lengths show saturation 
around one bond distance at low temperature, demonstrating that the system is 
in the quantum disorder regime at these three field values. 



V^Ej=i(5'|+j'5'|) and the magnetization (5) using the Trotter-Suzuki quantum 
Monte Carlo technique (QMC), where () is the thermal average and N is the total 
number of spins. We apply the cluster algorithm in the imaginary time direction, and 
the heat bath algorithm in the spatial directions. We warm up the system using 2 x 10^ 
Monte Carlo steps (MCS), and perform the measurement with another 3 x 10^ MCS. 
The thermalization is checked to be good. We take up to 64 ensemble averages using 
parallel computation. We use the periodic boundary condition in 3+1 (imaginary time) 
directions, and compute the spin correlation length along the z-direction, which can 
be summarized by the function ^(n, L, T, K), where n is the size in the imaginary time 
direction, L is the size in the x and y directions, and T is the temperature. There are 16 
spins in the unit cell of the pyrochlore lattice. Our QMC results have been calibrated 
and agreed with the exact diagonalization results of one unit cell. In order to obtain the 
correlation length in the thermodynamic limit ^(T, K) = lim„_j.oo, L-yoo ^{n, L, T, K), 
scaling analysis in the imaginary time and x and y direction is needed. To achieve 
this, we compute for lattice sizes LxLxz of 2x2x3, 3x3x3, 4x4x2, and 5x5x2 
unit cells, which contain 192, 432, 512, and 800 spins, respectively. To reach infinite 
n limit, we compute 4 different n, up to 40000 for L = 4, and up to 30000 for L = 5 
for the smallest T, and apply the finite-sized scaling analysis to obtain £^{L,T, K). 

In Fig. ([l^), we show the finite size scaling of the correlation length ^{L,T,K) 
as a function of 1/L. We compute for three different values of the transverse field, 
K = 0.05 J, 0.025 J, and 0.005 J. We only plot the results of the lowest temperatures 
that we are able to reach as the examples. They are as low as 2.5% of K for K — 0.05 J, 
and 5% of K for K = 0.025 J and 0.005J. In many regards, these are low enough to 
show the ground state properties, since K is the energy scale in the spectrum of 
low energy. All spin correlation lengths in Fig. ([T^) are below one bond distance 
from L = 2 to L = 5, and do not show much size dependence. First, we observe 
that the spin correlation in the z-direction shows exponential decay with correlation 
length which does not change with z. Second, the spin correlation shows little size 
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dependence in the x and y direction. Therefore, from the results of the short-range 
spin correlation and the scaling relation as shown in Fig. (jl^), we concluded that our 
calculations already represent the thermodynamic limit. ^(T, K) is computed from 
S,{L, T, K) by taking the mean of the values for i = 4 and L — 5. We plot ^(T, K) as 
a function of T for those field values in Fig. ([Tjs) . The correlation length shows neither 
significant temperature dependence nor K dependence, which is very similar to the 
kagome case [12l [13]. It remains very short ranged down to T equal to 2.5% of the 
field. Therefore, the ground state remains in the spin disorder phase. The numerical 
error in our calculations comes mainly from the ensemble average. The largest error is 
below ±3% down to very low temperature. The temperature dependence that shows 
the saturation of the spin correlation length implies that there is no phase transition 
down to zero temperature. As a first result, the system is a spin disorder for small 
field, i.e. disorder from disorder. 




Figure 2. (Color online) (a) A cartoon picture for the string excitation. Flipping 
spin at site ii generates a monopole (green ball) and anti-monopole (blue ball) 
pair. Further flipping at site 12 makes the anti-monopole hop. In this figure, the 
anti-monopole hops between 5 tetrahedra following the blue dash line, which is 
the flux-tube excitation of the energy ~ I — j- . The line in brown is an example 
of the pyrochlore hexagon which is the shortest path by which monopole and the 
anti-monopole can annihilate, (b) The phase diagram of the antiferromagnetic 
quantum Ising model in the pyrochlore lattice at zero longitudinal magnetic 
field. Deconfinement only occurs at K = 0. With non-vanishing K, the system 
becomes the quantum confined phase and undergoes a continuous transition to 
the paramagnetic phase. 



The next question is whether the classical monopole-anti-monopole pair is 
confined or deconfined in the presence of the quantum perturbation. Let us consider 
the following string operator [HI |3j 

Oi^ S+S'S+S- ■ ■ ■ S+ (2) 

where ik labels the sites of a connected string without intersection. As mentioned in 
the introduction, Oi excites the a monopole and an anti-monopole at each end with the 
separation distance I in the unit of bond length. An example of the string excitation 
is given in Fig. ([2^). They can be chosen in a random- walk way, but restricted to 
the condition that one tetrahedron does not contribute more than one bond to the 
string excitation. The string will never intersect but can be open or closed. A closed 
string corresponds to the annihilation of the monopole and anti-monopole, leaving 
a flux-tube-like excitation. In the classical Ising model, a flux-tube does not cost 
energy. In addition, we consider the string length I to be odd since [M^,Oi] = Oi, 
where = Suppose {Q\OjM^Oi\n)/{n\OlOi\n) = m^, where |f2) is the 
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ground state of finite K, the monopole and the anti-monopole sliarc ~ i.e. 
fractionalization. 

Now, let us consider the open string case and compute the energy of Oi by the 
following equation of motion 



-[[Oi,H],H] 

- J^UiOi - JKU2O1 
I 



' Sa • •■ S: 



fc=l 



fe=i 
1 ' 

-K^ s+s- ■■■{S+ - s-)-- s+ 



k=\ 



where 



27^2 ^ (- 1) Si S--- Sl^ ■ ■ SI- • • 5+ , 

m<n 

= Sl,^{l-\))\ 

=H5^+z5| ^St^, 

Us = S!,~il~l), 

U4k^U3 + {-l)''-^S^,+hk, 



(3) 

(4) 
(5) 
(6) 

where /i^ = 1 for fc = 1 and I while hk = 2 for others, and label the nearest neighbor 
spins of the ik site. In our system, one site has 6 neighbors. So far, the equation of 
motion as given through Eq.([3| to Eq.Q is exact. We compute the energy in the 
K <^ J limit by applying mean field theory. Take the following expectation values: 

{Ui)o:^{{-S!,~~S^,f)o,^Ao, 
{U2)o, - 0, 

" (-5f,-5f; + iSf. + i)o, fork = l ] 



(-S^,-S^, + ^Sf, + ^)o, fork = l 
(^S^^S^ + \{^\)'^St^V)o, 
for k = 2,3,..,/- 1, 



Pk, 



(8) 



where ()o^ denotes the expectation value of one string state, i'l labels the sites of the 
first tetrahedron without ii, and likewise for i[. Eq. ^ becomes 



df^ 



1. 



fe=i 



7nn.l-~2 



(9) 
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where i_i denotes the string operator of bond length I — 1 without 5"; in the Oi, 
and i_2 are the ones of bond length I — 2 without 5'^'^^ and S^^ in the Oi. 

E g. |9]) is noth ing but the simple harmonic equation for Oi with the eigenvalue 

uj = ^ J^Aq + ^ K"^ coupling to the strings of length / — 1 and 1 — 2. It can be solved by 

computing the equations of motion (EOM) for O' k,i-i and 0^„ ;_2' ^'^'^ obtain the 
set of EOM for the strings from / to the unit lengths. Solving the set of EOM bottom 
up, the right hand side of Eq. ^ becomes a large number of the quantum fluctuations 
of the string excitations of all lengths from Z — 1 to unity. For example, in Eq. Q 
there are I O'l^ and l{l — l)/2 O'J^^ ;_2 fluctuations, and they have the equations of 
motion of / — 1 string fluctuations of Z — 2 length, and (Z— — 2)/2 string fluctuations 
of / — 3 length, and so on. One can safely assume that their phases are random so 
that the inhomogeneous part in Eq. ^ is not crucial. Then, the string excitation O; 
is the quasi-eigenstate with the energy uj for small K. The procedure given above is 
the generalized random phase approximation for the string- like excitation O;. 

The monopole-anti-monopole pair is massive, because a single spin excitation 
takes the system away from the ground state manifold, which corresponds to the 
string excitation of unit length. The mass gap of the monopole-anti-monopole pair 
is Jy/~Ao — J ioT K ~ 0. Because there is macroscopic degeneracy in the ground 
state manifold, and also because the closed flux tube connects different ground states, 
the system is deconfined at the classical level. Our result of w is consistent with this 
scenario. For finite K, uj is /-dependent. It can be understood as the mass gap from 
the J term and the binding energy from the K term. When K is small, the binding 
energy is proportional to I ^ . As the binding energy is proportional to the separation 
of distance, the monopole-anti-pole pair is in a confined phase in the presence of the 
quantum perturbation. It is because the operator Oi actually is a projection operator 
for spin-1/2. As the ground state is spin disorder, namely < >= at every site, 
the operator creates a sequence of ordering sites which is nothing but the Dirac string 
of which the energy is proportional to its length. 

We note that the slope of the confining potential is small when K <C J. It 
is nearly deconfined. Although the open string excitation is an energetically costly 
quasi-eigenstate, the fiux-tube excitation of the closed string is less expensive. In the 
pyrochlore lattice, the shortest closed string is the pyrochlore hexagon the energy of 

which is y/3K. In general, the energy of the closed string is \J\k. Furthermore, 
we remark that the mechanism of the quantum confinement by the transverse field 
is non-trivial. As it is in the one-dimensional ferromagnetic Ising model, spinous are 
deconfined in the classical case and remain deconfined when the transverse field is 
turned on. It is the longitudinal field that makes it confined. 

For K ^ J, the mean field values of Eq. Q are no longer valid. The terms 
dominate in the EOM, and O/ is far from being the eigenstate. However, we can still 
compute the energy by {<^\0\HOi\'^) / {'^\0\0i\^) , where |$) is the ground state in 
the large K limit. It is not hard to obtain the leading term of the energy \IK. It is 
not surprising that the paramagnetic phase is a confined phase. Therefore, this mode 
is confined for any finite K, which implies that the paramagnetic phase in small K 
undergoes a continuous transition to the one in the large K limit. In other words, the 
phase in the small X is a cooperative paramagnet. 

Finally, we compute the effective Hamiltonian for small K using the Rayleigh- 
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Schrodinger degenerate perturbation theory, which is the following: 

where Jring = ^ 6 J ° ; Jfc labels the k*'^ site of the hexagon j, and the second summation 
is over all hexagons in the pyrochlore lattice. An example of a hexagon is given in 
Fig. ([2^). Eq. (10) is the classical Ising term plus a ring exchange term, which was 
previously studied by Hermele et al. [5]. They found the deconfined phase of the 
emergent pyrochlore photon in the large Jring limit. The present work fills the blank 
in the small Jring region. We note that the large Jring limit is different from the 
large K limit. Nevertheless, although Eq. (10 1 describes the physics of small- if, the 
open monopole pair is not in its spectrum, since it is massive in the order of J. For 
a close string, the Jring term is proportional to its rotational energy moving around 
the pyrochlore hexagon. 

In summary, our work introduces a new cooperative paramagnet in the three- 
dimensional pyrochlore system. The antiferromagnetic quantum Ising model in the 
pyrochlore lattice is a spin disorder ground state, and the monopole-anti-monopole 
pair is confined. The deconfinement only occurs at the K = point. Unlike the 
perturbation of the classical dipolar interaction, the quantum perturbation by the 
transverse field introduces the quantum confinement. The disorder phase adiabatically 
connects to the paramagnetic phase in the large field limit. Some interesting research 
directions include the finite temperature properties, the computation of experimentally 
measurable quantities, and search for new quantum spin liquid in other 3-dimensional 
systems. 
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